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' Abstract 

\ We examined the changes in swimming behaviour of the bacterium Rhodobac- 

O ' ter sphaeroides in response to stepwise changes in a nutrient (propionate), 

following the prestimulus motion, the initial response and the adaptation 
ff^ I to the sustained concentration of the chemical. This was carried out by 

' tethering motile cells by their flagella to glass slides and following the ro- 

O ■ tational behaviour of their cell bodies in response to the nutrient change. 

ir\ , Computerised motion analysis was used to analyse the behaviour. Dis- 

tributions of run and stop times were obtained from rotation data for 
tethered cells. Exponential and Weibull fits for these distributions, and 
variability in individual responses are discussed. 
^ ' In terms of parameters derived from the run and stop time distribu- 

■ tions, we compare the responses to stepwise changes in the nutrient con- 

centration and the long-term behaviour of 84 cells under twelve propionate 
concentration levels from 1 nM to 25 rriAI. 

We discuss traditional assumptions for the random walk approxima- 
tion to bacterial swimming and compare them with the observed R. 
sphaeroides motile behaviour. 



Keywords: bacterial chemotaxis, Rhodobacter sphaeroides, flagellar motor, 
Poisson process. 
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Introduction 

Bacteria populations exhibit intricate spatio-temporal patterns as a result of 
reacting actively to environmental cues, such as nutrient concentration, light 
levels, oxygen concentration, osmolarity, pH, presence of toxins or metabolites, 
electrochemical stimuli, etc. Individual cells integrate these signals to produce a 
balanced response, resulting in active locomotion. 

Bacterial locomotion in general can be achieved by different means, but most 
species swim using flagellum-driven motility. Rhodobacter sphaeroides is a purple 
non-sulphur photosynthetic bacterium. It has a single subpolar flagellum, which 
propels the cell forward, rotating only clockwise. When the flagellar motor stops, 
the functional helix relaxes to a short wavelength, large amplitude structure and 
Brownian motion appears to cause reorientation of the cell body as described in 
0. When the motor starts to rotate again, a functional helix reforms and the 
cell swims in a new direction. Bacterial motion usually can be described as a 
random walk with a good approximation. From the mathematical point of view 
it is often assumed that times of runs and stops are distributed exponentially. 

A wide range of phenomena occur when a signalling cue is non-uniformly 
distributed in space and/or time. Such stimuli can lead to directed or biased 
motion of cells, i.e. can cause a tactic response P]. This study is a step towards 
understanding the physical nature of taxis in bacteria subjected to changes in 
chemoattractant concentration. We start with a description of the experiment 
and analysis of data and then present results of fitting times of runs and stops 
with a generalisation of the exponential distribution. 



Description of the experimental setup 

We are exploiting data for temporal sensing in R. sphaeroides, which has been 
investigated using tethered cells and stepwise changes in the concentration of 
proprionate. The data were collected by H. Packer (Microbiology Unit, Depart- 
ment of Biochemistry, Oxford University). 

Growth conditions. Rhodobacter sphaeroides WS8N (a wild type nalidixic 
acid resistant strain) was grown at 30°C aerobically in the dark as described in 
(author?) [4]. 

Tethering. In order to analyse the behaviour of individual bacterial flagel- 
lar motors, cells were tethered by their fiagella using antifiagellar antibody in a 
flow chamber and the behaviour of individual motors was analysed using motion 
analysis. Cells were harvested by gentle centrifugation, washed and resuspended 
in nitrogen 10 mM Na HEPES buffer (pH 7.2) containing 50 g/ml of chloram- 
phenicol, aerated by shaking at 30°C. The cells were starved for 45 mins and 
tethered in the flow chamber as described in 0, |^ . 

Motion analysis. The tethered cells were viewed via a phase-contrast micro- 
scope with an attached video camera. Measurement of cell rotation was made 
using the AR0T7 software on a Bactracker. The software can analyse up to 10 
cells per field. Cells that were rotating without touching others were measured. 



A. A. Chernova et al: Response kinetics 



3 



Data points were taken at the video frame rate (50 Hz for interlaced images) 
and the raw data were downloaded as ASCII files for analysis. 

Tethered cells were subject to stepwise addition of propionate at the twelve 
levels with concentrations from 1 nM to 25 mM (i.e. from sub-saturating to 
saturating concentrations) over 3-5 minutes time interval. 

Data analysis 

During the experiment rotation speed of a single tethered cell was recorded every 
0.02 s. Typical time series of responses are shown in Figure ^ The raw ASCII 
data files were processed to eliminate system noise, each point in the curve was 
replaced by the weighted average of its nearest nine neighbours by the method 
of Savitzky-Golay described in [S]. An advantage of this filter is its sensitivity, so 
that subtle changes in bacterial motor dynamics are preserved in the data. The 
noise arising due to vibrational forces acting on a stopped cell must be treated 
separately. 

A typical distribution of rotation speed data (pre-stimulus period) is shown in 
Figure El The distribution is bimodal. The first peak around zero corresponds to 
stops of the motor; the negative values of the rotation speed appear when the cell, 
subject to vibrational forces, turns in the opposite direction. Vibrational forces 
contribute in a different manner to the distribution of speeds around the second 
maximum. A physical reason for this is the small value of the Reynold's number 
for a bacterium in liquid medium, which means that a moving cell is surrounded 
by a relatively more stable environment than a stopped one. The complexity of 
the frequency distribution of the rotation speed of a moving R. sphaeroides cell 
reflects the intricate behaviour of the bacterial motor. In comparison with the 
virtually constant-speed motor of Escherichia coli, in R. sphaeroides, the motor 
is known to be a variable-speed motor and the speed variation can be seen both 
in individual motors and between cells under the same external conditions. More 
detail on the stop/start mechanism of the motor in R. sphaeroides is given in 



For each cell the pre-stimulus records were used to define a threshold value of 
rotation speed. From the distributions of the rotation speed data we determined 
the two most frequent values (modes); the threshold ujq was chosen to be half 
of the maximal modal speed (Figure El). The cell was considered as stopped at 
time t if the detected rotation speed uj{t) was less than ujq, otherwise the cell 
was assumed to be rotating (running). 

Distributions of run and stop times 

Implementing the threshold criteria for the rotation speed data, we obtain se- 
quences of run and stop times for each cell. We computed mean values and 
standard deviations of run and stop times averaged over all the cells in one ex- 
periment (i.e. subjected to the same level of propionate) vs. the concentration 
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Figure 1: Experimental measurements of rotation speed of a single tethered R. sphaeroides 
cell. Vertical lines show the time intervals during which propionate was applied. Outside these 
regions the propionate concentration is zero. A delay is observed between the application of 
propionate and beginning of response. Duration of the delay differs from cell to cell. This part 
of data was not included in the analysis. Data from 1 /J.M, 100 fiM and 1 mM experiments. 
(Unpublished data). 



of propionate. Relative mean durations of runs and stops, defined as differences 
between mean duration of runs (stops) in the presence of the nutrient and mean 
duration of runs (stops) over pre-stimulus period, are shown in Figure El It was 
found that on the average the mean duration of runs remains approximately con- 
stant for cells which are subject to stepwise changes in propionate concentration 
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Figure 2: Determination of tlie thresliold level wq from the distribution of rotation speed 

data (pre- stimulus period). AH distributions of the rotation speed have two local maxima (two 
modes). The threshold level was chosen to be one half of the largest modal rotation speed. 
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Figure 3: Difference between mean values during application of propionate and pre-stimulus 
period vs. propionate concentration. 



from zero level to some rii, ni < 1 [iM. The stepwise changes from to some 
value ^2, 1 iiM < ^2 < 10 mM, cause an increase in the mean time of runs, 
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while as a result of the stepwise changes from to n^, > 10 mM, the mean 
duration of runs decreases. 

The mean duration of stops was not observed to change significantly in re- 
sponse to stepwise addition of propionate at the concentrations from 1 nM to 
500 nM. When the change was from to some concentration n4, n4 > 500 /iM, 
the average duration of stops increased, but the increase was not larger than 1 
s. 

Figure 0] illustrates how a single cell responds to stepwise addition and the 
following new constant level of the nutrient (here from zero (background) con- 
centration to 1 mM) : the cell tends to stop less frequently (the first mode disap- 
pears from the distribution of rotation speeds); duration of runs increases (up to 
4 times in this instance); there are less long stops than at the background level. 



Fitting with the Weibull distribution function 

In traditional modelling of bacterial motion it is often assumed that the tem- 
poral part of the random walk can be described by a Poisson process, or, more 
generally, alternating renewal process with exponential holding times, i.e. wait- 
ing times between runs are distributed exponentially (see, for example, 0,I3])- 
Here we attempt to justify this assumption by analysing distributions of both 
run and stop times obtained from the data. 

Distributions of run and stop times derived from the data were fitted with a 
Weibull distribution function, given by 

F^(t) = l-exp(-— ), (1) 
with probabihty density function 

fit). = -(-) exp(-^). (2) 

The Weibull distribution is a generalisation of the exponential distribution, in 
the sense that if /? = 1 then the former reverts to the latter with parameter 
ii = ri 

Fe(t) = l-exp(--). (3) 

A linear equation equivalent to (H} is 

$(r) = /3r - a, (4) 

with T = \nt, $ = ln(— ln(l — F)) and a = (3\nT]. An example of fitting data 
for a single cell is presented in Figure El The cumulative distribution function 
(cdf) calculated from data was converted into a form suitable for linear fitting 
with (jH). 
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Figure 4: a) Rotation speed (prestimulus period), b) frequency distribution of the rotation 
speed (prestimulus period), c) rotation speed (response to 1 mM propionate), d) frequency 
distribution of the rotation speed (response to 1 mM propionate). Percentage of e) runs 
(prestimulus period), f) stops (prestimulus period), h) runs (response to 1 mM propionate), 
g) stops (response to 1 mM propionate) as a function of time since the beginning of each 
run/stop. 

Results 

Figures IHl and [7| show the parameters rj and P estimated for each cell in all the 
experiments. 

As can be seen from the WeibuU parameter estimation, there is a significant 
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Figure 5: Fitting the Weibull function in the form to run and stop times of a single 
cell during pre-stimulus period (above) and after (below) addition of 1 mM propionate. 
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Figure 6: The Weibull parameter rj estimated for each cell in all experiments for run and stop 
times before and after addition of propionate. The upper plots show the parameters estimated 
in the absence of nutrient; the values on the horizontal axis correspond to different datasets, 
according to the subsequent application of propionate. For the three smallest concentration 
we used the same pre-stimulus segment of the dataset. Two values of rj are not shown in the 
plot rj = 65.04 estimated for run times at 25 mM and rj = 233.20 estimated for stop times at 
25 mM (a different cell). 



variability in the cell behaviour. 97.7 percent of all rj values lie within the range 
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Figure 7: The WeibuU parameter (3 estimated for each cell in all experiments for run and stop 
times before and after addition of propionate. When P = 1, the distribution is exponential 
with coefficient equal to the corresponding value of rj. The upper plots show the parameters 
estimated in the absence of nutrient; the values on the horizontal axis correspond to differ- 
ent datasets, according to the subsequent application of propionate. For the three smallest 
concentration we used the same pre-stimulus segment of the dataset. 



from to 12, and two values were found to be 65.04 (for run times, 25 mM 
propionate) and 223.20 (stop times, 25 mM propionate, different cell). 

The second parameter, (3, significantly deviates from unity for both run and 
stop times before and after addition of nutrient. 

One of the possible conclusions from analysis of norm of residuals of the fits 
is that the WeibuU distribution is a better description for stop times (both in 
the pre-stimulus period and after change in the chemical concentration) than for 
run times. 

The results of the analysis of run and stop times do not allow us to draw a 
single conclusion from the data, mainly due to a great diversity in individual cells' 
behaviour. Cases of successful fitting of a run/stop distribution with WeibuU 
functions are found as often clS CclSGS with poor results of such fitting. In some 
instances, like for the distribution of run times for a cell subjected to 1 mM 
propionate (Figure EJ, a linear combination of Weibull functions could be chosen 
for fitting. 

Preliminary analysis of independence of successive time intervals also show 
a high variability among individual cells, and does not allow us to talk about 
a "typical behaviour". On the average, however, the autocorrelation function 
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calculated for the sequences of run and stop times does not differ significantly 
from the autocorrelation function of a Poisson process. 

These complexities discovered in the data suggest that a detailed investiga- 
tion into the run and stop times distributions can produce some interesting 
results. We have chosen the WeibuU distribution function as a simplest gener- 
alisation of the exponential distribution. Such generalisation leads, in fact, to 
re-consideration of several modelling hypotheses for chemotaxis in bacteria. Ex- 
perimental motivation of the generalisation and modification of the temporal 
and spatio-temporal models is a subject of our current research. 
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